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I Abstract — The Casimir effect for a massless scalar field with Dirichlet and periodic boundary 

conditions (b.c.) on infinite parallel plates is revisited in the local quantum field theory (Iqft) 
framework introduced by B.Kay. The model displays a number of more realistic features than the 
^ ' ones he treated. In addition to local observables, as the energy density, we propose to consider 

d \ intensive variables, such as the energy per unit area e, as fundamental observables. Adopting 

this view, Iqft rejects Dirichlet (the same result may be proved for Neumann or mixed) b.c, 
and accepts periodic b.c: in the former case e diverges, in the latter it is finite, as is shown 
by an expression for the local energy density obtained from Iqft through the use of the Poisson 
summation formula. Another way to see this uses methods from the Euler summation formula: in 
the proof of regularization independence of the energy per unit area, a regularization-dependent 
surface term arises upon use of Dirichlet b.c, but not periodic b.c. For the conformally invariant 
scalar quantum field, this surface term is absent, due to the condition of zero trace of the 
energy momentum tensor, as remarked by B.De Witt. The latter property does not hold in 
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the application to the dark energy problem in Cosmology, in which we argue that periodic b.c. 
might play a distinguished role. 
0. Introduction and Summary 

The Casimir pressure of the electromagnetic field enclosed by (infinitely thin) parallel plates, 
measured by Spaarnay, is one of the most famous objects in quantum field theory (and in 
quantum optics) [I], [2]. A good up-to-date account, including recent experiments, is to be 
found in the volume Seminaire Poincare, vol.1, with the title "Energie du Vide". 

In spite of the well-known exact solution for the pressure (pp, [2]), the energy per unit area 
appeared to remain divergent, due to the (nonintegrable) divergence of the energy density at 
the boundaries - a phenomenon analyzed quite generally in the pioneering paper of Deutsch and 
Candelas [3] - until recent work by H. Kiihn clarified the situation, showing that the divergences 
due to the electric- field and magnetic- field components exactly cancel [Ij. This argument does 
not, however, hold for other fields, which may play a role in Cosmology, particulary in the 
problem of dark energy ([5], [6], [7]). There, the energy density - the (time)-OO component 
of the energy- momentum tensor Tqq{x) - is as important an observable as the pressure. If 
it diverges, or is ill-defined, as in the case of general fields enclosed by parallel plates, the 
situation remains highly unsatisfactory from a conceptual point of view. In fact, this is the most 
elementary example of cutoff-dependence of (in principle) observable quantities, which has been 
emphasized by Hagen in more general situations [8]. 

In this paper we revisit the Casimir effect for a massless scalar field with Dirichlet and 
periodic boundary conditions (b.c.) on (infinitely thin) parallel plates. Other b.c. (Neumann or 
mixed) may also be handled by the same methods, and yield results qualitatively similar to the 
Dirichlet case. 

In section 1 we introduce the general framework and ideas, which go back to B.S. Kay [9] 
and L. Manzoni, G. Scharf and one of us (W. Wreszinski) ([lOj. |llj). (See also the paper by 
Hollands and Wald for a very stimulating introduction to the local quantum field theoretic (Iqft) 
aspects of the problem [IS]-) There we show that Iqft yields a definite formula (1.27) for the 
energy-density operator. 

In section 2, we prove that closed form expressions (2.12a,b) for the true vacuum's ( i.e., 
in the presence of the plates) energy density follow from (1.27). This is the answer provided 
by Iqft to the problem. For Dirichlet b.c. this energy density (2.12a) displays, however, a 
nonintegrable divergence at the plates, thereby yielding (nonphysical) infinite values for two 
fundamental observables, the energy per unit area and the pressure. The renormalized one- 
point function is, moreover, also given by a nonintegrable function and , thus, does not define 
a (Schwartz) distribution. In contrast, for periodic b.c, the energy per unit area is finite and 
homogeneous by (2.12b). The reason for the infinity in the Dirichlet case is that only for 
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two space-time dimensions is the renormalization assumption (1.5) powerful enough to yield a 
divergence-free theory [9]. Otherwise, additional divergences arise, due to the sharpness of the 
surface, as remarked in the pioneering paper of Deutsch and Candelas [3]. The problem occurs 
whenever the attempt is made to impose b.c. on quantum fields (called "unnatural acts" by R. 
L. Jaffe in [IE]), i-e., to restrict quantum fields to sharp surfaces. 

In a more specific Iqft context, a similar problem arises in the restriction to a causal surface 
(the horizon) in connection with the problem of localization entropy [19] . The effect observed 
in [19] is, as remarked there, the thermal counterpart of Heisenberg's old observation that, in 
approaching conserved global charges as spatial limits of integrated charge densities, one must 
control the vacuum fluctuation in the neighbourhood of the boundary. The divergence of the 
inverse "split distance" in [19] has, as also remarked there, nothing to do with an ultraviolet diver- 
gence. Both features have strong analogy to the results in section 2. There, the ultraviolet cutoff 
drops out from the calculation: it is merely due to the introduction of the (cutoff-dependent) 
regularized flelds (1.22) as an intermediate step, which seems to be unavoidable due to the subtle 
cancellations occurring in (1.27), and was also the route followed by B.Kay in [9], who added 
that "experience in quantum fleld theory leads us to expect that, after renormalization, results 
will be regularization-independent for a wide class of regularization procedures." 

In this paper we are able to develop this conjecture, thereby improving on [9J in this aspect. 
This turns out to be a nontrivial task, but it provides an additional, important information: we 
prove in section 3 that the above-mentioned regularization independence (RI) within a wide class 
implies - in a purely mathematical context, of a rigorous asymptotic analysis along the lines of 
the Euler-Maclaurin summation formula [23] - the existence of a regularization-dependent 
surface term (RDST) (the second term in (3.17a)) which diverges as the cutoff is removed. 
This very nonlocal term is, as remarked, due to sharpness of the (infinitely thin) surface in the 
case of Dirichlet b.c. 

The RDST has a long history, both in the physical and in the mathematical literature. In the 
physical literature, it is associated with the names of Symanzik [IB], Candelas [33], Deutsch and 
Candelas [3], and Barton [34]. in section 4 we call them, for brevity, Symanzik counterterms. in 
the mathematical literature, the surface term is well-known as the area term in Weyl's formula, 
see, e.g., [35]. However, this is not intrinsic to Iqft, but depends on the regularizer (which is the 
exponential one (2.1) in [35]). The introduction of a regularizer as an intermediate step seems to 
be unavoidable because of the subtle cancellations occurring in (1.27). The regularized version 
of (1.27) becomes (1.28), but, as seen, the regularizer drops off in section 2. 

On the other hand, Elizalde [21] has shown that Hadamard regularization of the integrals 
occurring in the models treated in [18] yields terms which, in the limit of a sharp surface, remain 
finite. In section 4 we show that the same happens in the present model. This is accomplished 
by the use of the Poisson summation formula, whose importance is even better appreciated 
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in the massive case, where logarithmic divergences occur [36]. The interest of the Hadamard 
regularization is two-fold: it is, as Elizalde remarks, able to separate and identify the singularities 
as physically meaningful cutoffs associated to a finite width (see (4.6)), on the other hand the 
corresponding one-point function (which is a constant for periodic b.c.) is a pseudo-function 
with the same type of singularity with respect to the singular surfaces as the (old) Hadamard 
form (5.3) of the two point function. The latter has a fundamental meaning in locally covariant 
quantum field theory [31], and it seems of interest to inquire upon the eventual significance of 
this structure in the present context, see section 5, where we also present our conclusions. 

We end this introduction with some important remarks concerning the physical motivation 
of the present paper. The RDST are identical to those surface terms mentioned by B.S. de 
Witt in his elegant analysis [37] (p. 307) (see also L.H.Ford [39]). As he remarks, similarly to the 
electromagnetic case, for a conformally invariant massless scalar field, with traceless stress 
tensor, these surface terms vanish! The question may thus well be asked whether this is not a 
better and more fundamental proposal to get rid of the bothersome (divergent) surface terms, 
why is the scalar field of any interest at all. Perhaps the infinity just reflects the fact that the 
massless scalar field is an unphysical object, at least in the context of the Casimir effect. 

We now argue that this is not so, at least as regards the possible relevance of the Casimir 
effect to the dark energy problem in Cosmology, for which there is both theoretical [7J and 
experimental [S] evidence. This is our main physical motivation; see jT5] for an introduction to 
the subject. Of course, the geometry of the present model is not relevant to these (potential) 
applications: an adequate choice would be a closed surface involving a point, e.g., a cube. The 
parallel plates offer, however, a much simpler case in which the Casimir problem for a scalar 
field may be studied rigorously. 

The first point is that the stress tensor is not supposed to be traceless, being of the form 

T^,u = gtMuP (0.1) 

where g^u = (1,-1,— 1,-1) for flat (Minkowski) space-time - a good approximation for the 
present cosmic time |38j. Above, p denotes the energy density. We now come to the second 
(related) point. From (0.1), 

P = -P (0.2) 

where p is the pressure. Taking for p the dark energy density, and p the pressure exerted by dark 
energy, (0.2) seems to be very well satisfied by the data for dark energy [7]. For a conformally 
invariant massless scalar field, 

p<0 (0.3) 

|37j . Assumption (0.3) would contradict, together with (0.2) the fact that the observed pressure 
is negative, which is, perhaps, the most surprising property of dark energy, and the one which 
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led Turner [7] to declare that "vacuum energy is almost the perfect candidate for dark energy" . 
(The argument of [46], leading to the result that the energy density is nonpositive in quantum 
field theories, is not (generally) valid in the Casimir effect, because in [46] covariance of the 
fields with respect to the Euclidean group of translations is used in the proof.) It is, in fact, 
precisely the negative pressure that is difficult to obtain in alternative proposals for dark energy 
(see also [40]). (0.3) holds for the model treated in [9j. By (3.17a), it also holds for the finite 
part of the energy density of the massless scalar field, but the latter is positive for the massive 
scalar field for suitable values of the mass in the case of parallel plates [36], as well as for the 
inner problem of the cube p]. As an example of a massive scalar field , one may consider the 
hypothetical axion, which may have very small nonzero mass [41]. The latter is a candidate 
for cold dark matter [7j, and one may investigate its corresponding vacuum field. Of course, a 
massive scalar field has no conformally invariant counterpart, and the surface terms do occur 
in this case (except for periodic b.c, see section 5). We postpone further comments on this 
and related issues to the conclusion in section 5. We refer to pQ, [2] and [21j for (part of) the 
immense literature on the Casimir effect. For the applications to Gravitation and Cosmology, 
there is an early beautiful review by L. Parker [45], as well as recent papers by Milton, Elizalde 
and coworkers [42], [43], [44], where many other references can be found. 

1. General Framework 

We consider a massless scalar field ^>(x) on Minkowski space time x = {xq,x). The corre- 
sponding Hamiltonian density h{x) is given by (^ = c = 1): 



We also wish to consider the free (massless scalar) field restricted to the region K/^ between 
two (infinitely thin) parallel plates at z = and z = d: 




(1.1) 



Kd = <x = (xu, z), with xu = (x, y) £ and < z < d 



(1.2) 



with Dirichlet or periodic boundary conditions (b.c.) on the boundary dK^ of K^,: 




(1.3) 



The density operator corresponding to (1.1) is given by 
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H{x) =: h{x) 



(1.4) 



where the dots indicate normal (or Wick) ordering. Measuring H{x) is a local operation which 
involves only a small neighborhood N{x) of the space-time point x. Since, however, the state S 
of the system on is different from the vacuum state uj of (infinite) space-time even restricted 
to N{x) (see [9], appendix, for a discussion), the question arises: with respect to which state is 
the normal ordering (1.4)? In [9j, the following renormalization condition was imposed: 

io{H{x)) = (1.5) 



for all X in Minkowski space-time. This condition means that double dots refer to the infinite- 
space Minkowski vacuum state uj, and was motivated in [lOj . [llj by the fact that real boundaries 
consist of electrons and ions, and the field which interacts with them is quantized in infinite space, 
but one may also view (1.5) as an independent renormalization condition, as done by B. S. Kay 
in [S]. The assumptions of local quantum theory [T3] yield now a rigorous formula for S{H{x)) 
(see, again, the appendix of [9]): 

S{H{x)) = lim I- (dr,o,dxo2 + ds^dj;2)uJi^+ixi,X2) (1.6a) 

where 

CJ1,+ (X1, X2) = S($(X1)$(X2)) - C^($(X1)$(X2)) (1.6b) 



The scalar field of zero mass, quantized in infinite space in p dimensions, may be formally 
written 



$(p)(a 



1 



dPk 



a{k + a+{k)e'''-'= = $L^^(x) + $f (x) (1.7) 



where x = (xq, x),k ■ x = k^XQ — k ■ x, 



ko = LOj:= \k\ 



(1- 



and <!>_, <!>_(_ refer to the negative and positive-frequency parts in (1.7), i.e., those associated to 
a (resp. a"^), and satisfy 

(1-9) 



1 
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with 



1 fdPk 



(27r)P J 2uj 



-ik-x 



(1.10) 



Up, Up are annihilation and creation operators defined on symmetric Fock space over the (one- 
particle) Hilbert space Tip = L^(M*'), J^giTip), with a(/) antilinear, a~^{g) linear, such that 



;i.lla) 



on a dense domain of finite-particle vectors (see, e.g., [15], p. 208 ff.), where {f,g)'Hp denotes 
the scalar product on Tip. The vacuum Op is such that 



apimp = V/ G Hp 



(1.11b) 



The scalar field of zero mass on the region is formally given by 



(2vr)2 



[1.12a) 



where 



Ur {x)=e"'r-\\Un{z,D) 



(1.12bD) 



with 



Uniz,D) 



2 fn-K . 

-sin(-.) n=l,2. 



(1.12cd) 



and k\\ = (kr,, ky), x\\ = {x,y), 



ijj 



k«,n 



\h\' + 



/nvr 
W 



1/2 



for Dirichlet b.c, or 



Un{z,P) = \j\ exp(^^z) 



(1.12dD) 



(1.12bp) 



with 
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a;,- = 

kn,n 



2 /2n7r\2 



Kir) 



1/2 



n = ... ,-2,-1, 0,1,2,... 



(1.12dp) 



for periodic b.c. Above, Uj^^^ are (improper) eigenfunctions of (— A)^/^, where A denotes the 
Laplacean: 



(1.12e) 



The a, a"*" in (1.12a) are operator-valued distributions on J^si'H), where 



(1.13a) 



and 



oo oo 
= {/ I f{z) = l.i.m. CnUniz), with ^ |Cn|2 < Oo} (1.13b) 



n=l 



n=l 



where l.i.m. denotes the hmit in the topology of L^(0, d), and such that 

[a(/),a+(5)] ={f,9)nl 



(1.14) 



on J^Q, where {f,g)n denotes the scalar product on H- The vacuum is defined by 

a{f)nK, = \ff e H (1.15) 
Analogous definitions apply in the case of periodic b.c. The field has the two-point function 

(1.16) 



<^K,,_{X),^K,,+ {X )J = -D'j[J{x,X ) 



1 



where -D^^ is the distribution 



^k] {x,x') = (Xo - Xg, X, x 



dkn 



n=l fell 



.e^^l-(-l--;)[/„(2)C/„(z') 



(1.17) 



Due to (1.16), (1.17), the canonical commutation relations (CCR) are altered with respect 
to their free field values \i f,g ^ 5(R^) ^ C^(0, d) ® 5(M) (corresponding to the variables xy, z 
and xo), then 







(1.18) 
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not only when the points x = (xy, z, xq), x' = {x^^,z' ,x'q) in the supports of / and g are space- 
Uke to one another, but also whenever x is space-Uke to its mirror image (xj|,2(i — z' ,x'q), i.e., 
(x|| — + (z + 2; — 2(i)^ — (xo — Xq)^ > 0, corresponding to connecting the points x and x 
both by a ray and by one which suffers a reflection on the right plate (there is, of course, an 
infinity of other possibilities involving the left plate and multiple reflections). This is valid only 
for Dirichlet boundary conditions. 

See also |16] for the commutation relations for the electromagnetic field in the Coulomb 
gauge and Dirichlet b.c, and for a general discussion of boundary conditions in quantum 
field theory. 

We need also regularized fields. Let C : M — > M be a smooth function (1.19a) satisfying 

C(0) = 1 (1.19b) 
C'^''\x) dx < 00 VA; = 1,2, ... (1.19c) 

lim C('=)(x) =0 VA: = 1,2,... (1.19d) 



We introduce a class of regularizing functions C\ depending on a cutoff A with dimensions 
of length by 

Ca(^) = C{Auj^) (1.20) 

where lo^ is the frequency (1.8) (which in the case of ^k^, given by (1.12a), is given by (1.12d)). 
By (1.19a) and (1.20), 

lim CA(k) = 1 (1.21) 
The regularized fields and ^k^A corresponding to the formal and defined 

by 



^'^\x) = {27:)-^'' j-j=\ap{k)e-''-- + a+(fc)e*'=- 1 • Ca(A:) (1.22) 



00 



/V — > 1 f -* —iujf xc] 



— * — ZCJ-* Xn \ — # 

+ a+{k,n)Ur (x)e 'H'" • Ca(A:) (1.23) 
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where Ca is given by (1.19), (1.20), and has the property (1.21), as well as a normalization 
condition inherited from (1.19b). 

For the free field, we may, equivalently for (1.6), go from infinite space to a geometry with 
boundaries, by expressing the normal ordering (1.4) in configuration space by the point splitting 
technique which yields 

: = lim : ■ — 




*x dxQ dxQ 

" dxo dxo ^ i9xr°''^ ' ^ ' 

in the sense of operator-valued tempered distributions, where, in (1.24), 

$ = (1.25) 

Finally, from (1.1) and (1.3), 



H(x) = lim < - 



d^(x)d^(y) , ^ 

dxQ dyo 



+ 



which we take to be the Hamiltonian density (operator) describing the field, both free and with 
boundaries, in agreement with (1.5) and the discussion following it. In case we wish to describe 
the field with boundaries, the first three terms in (1.26) must be defined on symmetric Fock 
space on the adequate (one-particle) Hilbert space TC, i.e., the concrete representation of the 
field operator is dictated by the geometry. In the case of parallel plates, i.e., with the fields 
defined on K^, given by (1.2), with Dirichlet b.c. on dK^, given by (1.3), the field is (formally) 
given by (1.12), on J^sCH), given by (1.13). 

Thus, (1.26) has to be an operator on J^si'H), and therefore must be represented in the form: 



where D^x^ is given by (1.17) and the semicolons in (1.27) denote normal ordering, with respect 
to the emission and absorption operators in (1.12a) (which satisfy (1.14)) and the vacuum ^k^j 
defined by (1.15). 
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Taking into account (1.22), (1.23), we define the quantities corresponding to (1.1), (1.10) {inp 
dimensions), (1.27) and (1.17), denoting them by h^^\x), ^q'a ' ^A,Ka{x) and D^k] ^ respectively. 
H^ Kd represents the regularized Hamiltonian density in a theory with boundary dK^i given by 
(1.3). 

Let, now, p = 2, and Tir (resp. Tii) denote two copies of L'^{M?) corresponding to the right 
(resp. left) plates (components oidK^), J^s{T~(-r) and J^sCHi) the related symmetric Fock spaces, 
with vacua Qr, ^s, satisfying (1.11b) (with p = 2). The related Hamiltonian densities will be 
denoted by h)^\.{x) and h^\{x). They will make their appearance only in section 4. 

Further, let 5r \ denote the delta-functions associated to the right and left plates (see 
|24j . Ch.3, § 1, for delta-distributions and other singular functions associated to a regular sur- 
face) . 

Finally, let E^'^ denote the vacuum energy density corresponding to H^^k^- By (1.27), it is 
(with defined by (1.15)): 



Er{x) = {nK„HA,KA^)nK,) 

= - hm i ^ 

I y^x I (JXq 



(1.28) 



2. Poisson's summation formula and the Casimir energy density 

In this section we use for convenience a special regularizer, 

C{x) = e""^ x>0,xGM (2.1) 

which satisfies the conditions (1.19a), (1.19b) and (1.19c). In the next section we shall 
demonstrate that there is no loss in choosing any regularizer defining the above conditions for 
the calculation of the Casimir energy density. Putting (2.1) into (1.28) we obtain 

El--{S) = -^IA^ L3A;e~^['"""'-^"^"^^^]^° xe-^'=o- 



2 M 1 (27r) 



jdki jdk2 Y,\Un{x)\^ X e-^'^A (2.2) 



with, for the parallel plates, Un{x) either of 



UnM^) = ^\Jl sin(^z)e^('=^^+'=-^)n = l,2,... (2.3d) 



or 

1 / 1 / 2imT 



UnA^) =7f\/^exp(^z)e^('=-"+*^-^)n = ... ,-2,-1,0,1,2,... (2.3p) 



11 



where uj'^ = kf + k'^ + (if )^ Dirichlet b.c. and o;^ = fcf + + (^)^ for periodic b.c, 
using now and henceforth the letters D and P to denote Dirichlet and periodic b.c. We first 
consider Dirichlet b.c. From (2.2) we may write 

where k = \k\, and 

m^ji^lpjdki jdk,s^n\^z)e-'-- = 

^ +00 „ „ 
= E y^^i Jdk,s^n^^z)e-^-- (2.5) 

The summation term in the energy density in (2.4) can be extended to accommodate negative 
terms in n, noticing that the parity of integrand is even and the term n = yields zero. Let 
now, by Poisson's summation formula (PSF) (see, e.g., [20]), 

+00 +00 

^ /(2vrn) = fi^) (2-6a) 



where 



and 



/2tt 

n=— 00 ' m=—oo 



f{m) = ^= / dke-'^'f{k). (2.6b) 
V 27r 



/(27rn) = sin^ C^z) e"^"^" (2.6c) 



we obtain from (2.5) 



2 dK 2(27 



X W(|z)e-^V'^^^'^^n-; (2.7) 



Performing a change of variable k^ = k^/2d on (2.7) we get 

I = h + h (2.8) 

with 

1 +00 /•+CO p+00 r+oo I 

h = dki dk2 / dk3 ^kf + A;2 + fc2g-A^/fc?+fcI+fcig-i2m*3 (2.9) 

\ / ,™ — — OO "/ — CO — OO 
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2(27r) 



^ -i-oo /■+00 n+oo r+oo , 

Yl dki dk2 / dks kj + kj + k. 



m=— oo 



(2.10) 



By (2.4), (2.8), (2.9) and (2.10) we obtain 



+ 



(2vr) 



1 J 4 (2z)2 - 3A2 

2 I (2^)2 [A2 + (2z)2]3 

[2(m(i + z)] ^ - 3A2 
[A2 + [2(md + z)]2]^ 



(2m(i)2 - 3A2 
(2^)2 [A2 + (2m(i)2]3 



E 



+ 



+ 



[2(md 



3A2 



[A2 + [2(md - z)f 



(2.11) 



where the first term above is the term with m = in /2- The m = term in 1\ cancels exactly 
the first integral in (2.4). By (2.11), 



1^1 OO 

'2(2^ E (^^ ^ 4(2^ E 

^ ' m=l ^ ' ^ ' m=l 



1 



+ 



1 



(2.12a) 



(md+z)^ [md — z)'^ 

The first term in (2.12a) diverges on the plate where 2 = 0, the last term (for m = 1) diverges 
on the other plate where z = d. These dominant terms can be identified as stationary points in 
the PSF by (2.7): 

Indeed, writing sin?(^^z^ = 2i — ^) ™ (2-7), these points correspond to [(md ± 

z)li^ = 0, which lead to z = (for m = 0) and z = d (for m = ±1). 

A general method combining the PSF with the stationary phase method for obtaining asymp- 
totic estimates was developed in [20] (see also [22]). 

We now consider periodic b.c. in this case it is readily seen from (2.3p) that /(A) in (2.5) 
is given by the same expression on the r.h.s. of (2.5) with, however, sin(^z)^ replaced by one. 
This implies that I2) = in (2.8), and Ii is given by the same expression on the r.h.s. of (2.9), 
but with exp(— 2im(iA;3) replaced by exp(— imcZfcs). The final result is 



E'""'{z,d;P) = lim ET''{z,d;P) 



1 



rC(4) 



TT 



Thus, (2. 12) is the result provided by Iqft for the original problem. We define 

E^'^'iz, d) = for z <0 01 z > d 



(2.12b) 



(2.13) 



corresponding to the fact that the outer Casimir problem yields zero contribution to the energy 
density and pressure ( Remark 3.7 ). 
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Consider, now, a compact subset of (defined by (1-2)): 
i^f ^(f||,z),withf||(x,y)e [-:^,:^] X [-:^,^] andO<z<d| (2.14) 

with 

L1L2 = A (2.15) 



Let {xl^l (•)}„=! 2 be a sequence of smooth functions approaching, as n — > 00, the char- 
acteristic function of K^, and such that 

x5^1(x) = 1 for f e K^^'O U dK^ Vn=l,2,... (2.16) 

d 

where K'^'^ denotes the interior of Kf. 

We define the vacuum energy E'"°''^{A) of the region by 

hm hm [d'xx^;^\{x)E^,Z^^{xo,x)^E^-'{A,d) (2.17) 

whenever the double-limit on the l.h.s. of (2.17) exists, is independent of the regularization 
(1.19), (1.20) and does not depend on xq. The vacuum energy per unit area e - a basic observ- 
able quantity ~ is, then, defined by 

E'"<^^{A,d) 

e = lim 2.18 

A^oo A 

whenever the hmit on the r.h.s. of (2.18) exists. By (2.12a), E'""''^ (z , d; D) is not integrable over 
< z < d, and therefore 

ej^ = -00 (2.19a) 

However, by (2.12b), for periodic b.c, 

sp = -r (2.19b) 

90^3 ^ ^ 

The above differs from the finite part (3.17b) of ed by a factor 2~'*, due to the mentioned 
replacement of the summand 2m in exp{—2imdk3) by exp^—lmdk^) in (2.9). 

(2.19a) means that Iqft rejects Dirichlet (Neumann, mixed) b.c. in the Casimir effect for 
a scalar field (see also the last paragraph of section 5). For the moment, we just remark that 
this is consistent with the fact that formula (1.6) has being arrived at by the assumptions of 
local quantum theory ( [9] , appendix) , in particular the assumption that the states S and co are 
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locally quasi-equivalent (|14j. pg.l31). This is not true everywhere, there are singular points or 
surfaces (in our case in R^) in the neighborhood of which even such local quasi-equivalence must 
fail to hold, otherwise the renormalized one-point function (1.6b) uji^^(x,x) = uj^^i(z) which 
depends only on z (the two equal time variables being arbitrary), would exist as a Schwartz 
distribution by the ingenious discussion in the Appendix of [9j, which we have seen not to be 
true by (2.12a) in the case of Dirichlet b.c. (a similar proof holds for Neumann or mixed b.c). 
Local quasiequivalence is essential to thermodynamics and of great importance in the context of 
quantum field theory in curved space-time (which means that it should also play an important 
role in cosmology), see the discussion and references in |14j. 

In the next section 3 we provide an alternate picture of why Dirichlet and periodic b.c. yield 
results which are so qualitatively different in the context of cutoff (regularization) independence 
of e, where we show the emergence of the RDST mentioned in the introduction for Dirichlet but 
not periodic b.c. (remark 3.6). Notice that the forthcoming (3.17a,b) are consistent with (2.19a)! 
Section 4 relates these surface terms to Hadamard regularization of (2.12a), in agreement with 
previous work on the subject, also mentioned in the introduction. We leave to the last three 
paragraphs of section 5 (conclusion) a brief but important discussion of why we believe that the 
message of Iqft, our starting point, is clear and correct also in the context of the Casimir effect. 
3. The cutoff independence of the energy per unit area 

We now write (1.28), with a general cutoff satisfying (1.19) and (1.20), as 
A 



2(27r) 



2 



poo 

2d / dk k^C{Kk) + 



2.t[ikk,ir4)\k^-c{^,ir4^+^^]\ (3.1) 



d J \ \l \ d 



Rewriting the second term of the r.h.s. in (3.1) with a change of variable we have 

d 



{2ny 



poo -I ''I 

/ dk k^C(Ak) + — lim y g(m) } (3.2) 

m=l J 



with 



g{m) 



du^/uC{K^/^) (3.3) 

(rnvr/d)^ 

Now, we can introduce the Euler-Maclaurin sum formula, which yields ([23], p. 326) under 
assumptions (1.19c) and (1.19d): 
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^ 5(m) dq q'C{Aq) - -g{0) ^ 



m=l 



where 



with 



and 



1 f°° 

Sfc = -5fc(0) - j^j-^^ dti;2k+2it)9^'''+'Ht), A: = 1,2,... 



= (/'fc(i) mod 1 

where (pk can be obtained as follows 



X 

p.^ - 1 

n=l 



Theorem 3.1 



Under assumptions (1.19c) and (1.19d), (with /c = 1, 2, . . . ) is independent of /c, 

Sfc = Sfc+i V A: = 1,2, . . . 

Proof . By (3.5) 

= -5fc+i(o) - ^^^^ dtv'2fc+4(t)9(''+')(t) 

The first term on the r.h.s. of (3.10) may be written 

k 

^ (20! ^ ^ ^ [2(/fc + l)]!^ ^ ' 

By (3.7) and (3.8) it follows that (see [23], pp.320-321) 

MO) = 0, 
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and 



(3.12b) 
(3.12c) 



Rewriting the second term on the r.h.s. in (3.5) using (3.12c) with m = k + 1 and integrating 
by parts using (1.19d), we have: 



Sfe+i + 5^+1(0) = 



(2A; + 2)! 



1 /-CX) 

+ (-l)'=S,+i5(2fe+i)(o^ 



-'k ■ 



(3.13) 



Using (3.12b) above, with m = k + 2, and integrating by parts again we have 
which yields 

Finally, putting (3.15) and (3.11) into (3.10), we find 

S,+i = -5,(0) - ^^^^ dt^2k+2{t)g^''-''^ = 



(3.16) 



Theorem 3.2 

Let in addition to (1.19c) and (l.lQd), the normalization condition (1.19b) hold. Then 

rpvac 1 rca 

^ = J duu^C{u) +£/ + 0(A2) (3.17a) 



A 

where 
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with Ef independent of the regularizer, if the latter satisfies (1.19). 
Proof . By (3.2), (3.3), (3.4) and (3.16): 



jpvac 

A 



1 °° 1 
hm — y g{fn) = — 



m=l 



2^(0) + S2 



(3.17c) 



where, by (3.5), 

^2 = -2 ^ V-/ . 24^ 6! 



''^^("^(O) + - i r Mt)g('Ht)dt (3.18) 



By (3.3), 



1 1 1 



^(0) = - I du v/^C(A0I) = ^ I du u'^C{u) (3.19) 
2 2 7o A'^ 7o 



Again, by (3.3), 

-g^'Hm) = -2(— ) C[A—j (3.21a) 



By (3.18), (3.19), (3.21) and (1.19b), 

= -f g)%0(A^) (3.22) 

Putting (3.19) and (3.22) into (3.17c) we obtain , using B2 = (3.17a), with Sf given by 
(3.17b). 

Remark 3.1: It is remarkable that in (3.4) the second term in the l.h.s. is the vacuum term, 
and the third one, the surface term, appearing in a natural way as necessary subtractions in a 
purely mathematical context. This surface term will be reconsidered in the next section. 

Remark 3.2 : Usually the result is presented informally without the important last term in 
(3.5), and assuming that C satisfies C^''\0) = for all A; > 1 besides (1.19b), which is not 
satisfied by the special and important choice (2.1) ([27|, p. 138). See, however, ref. [28] for a 
much nicer approach to the subject. 

Remark 3.3 : S/j {k > 1) is called the (9^, 0) sum of the (divergent) series ^^^=1 9ij^)i where 
9^ refers to Ramanujan and to the reference point (the origin in our case). 

Remark 3.4 : Other classes of regularizers do not necessarily lead to regularization indepen- 
dence. See [9] for thorough discussion of several types of regularization. We believe the present 
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class is "natural" from the physical point of view, because it simulates a dielectric constant with 
suitable behavior at high energies. 

Remark 3.5 : The idea of the proof of Theorem 4.2 appeared in the second reference under 
Theorem p. 319, but only an incomplete sketch was given there. More importantly, however, 
the term ^^(0) was (wrongly) asserted there to contribute only to the A-independent terms, while 
it is precisely this term that yields the surface contribution. 

Remark 3.6 : The surface term in (3.17a) is absent for periodic b.c, because the latter allows 
for the m = term in (3.17c), which cancels it exactly. This explains the result of [lOj, which 
bears some similarity with the also very special model in |9j. 

Remark 3.7 : The external Casimir energy is zero, see |TTj, whose proof remains unaltered. 
Result (3.17a) is one-half of the result for the electromagnetic field, due to summation over the 
two polarization states in the latter. Polarization does, however, play a major role in explaining 
the cancelation occurring in [4]. 

Remark 3.8 : The series (3.17a) is a divergent asymptotic series, but the rest in (2.13) is - 
by (3.17c) - bounded by constant A^, so that its limit as A ^ 0+ exists and is zero. 

Remark 3.9 : The present approach also works for the inner problem for the cube, see [6] 



4. Symanzik counterterms versus Hadamard regular izat ion 

In the previous section we saw that RI of the energy per unit area is not true in the Dirichlet 
case unless a regularization-dependent surface term is subtracted from the (regularized) en- 
ergy density: the second term on the r.h.s. of (3.17a)! This term corresponds exactly to the 
replacement of the regularized vacuum energy density (2.2) by 



^ren,A 



1 



(2) 



(x) 



(4.1) 



The additional terms in (4.1) may be called "surface renormalization counterterms". they have 
been introduced by Symanzik [13j in a different framework of quantum field theory in the 

(2) 

Schroedinger picture and including interactions. In (4.1), the Hamiltonian densities h)^ j.{x) and 

(2) ' 

h)y i{x) and the vacua fir, have been defined in section 1, and correspond to the Hamiltonian 
densities and vacua of regularized quantum fields in (infinite) two-dimensional space (p = 2 in 
(1.22)). Thus, in (4.1), 

1 27r 



(f2j. , /i^ ^ J^r) 
and similarly for the left index. Finally, 



2 (27r)2 



dk k^C{Ak) 



rpvac 

ren,A 



hm ld^XX^^\{S)ErenA 
n— >oo J 



(4.2) 



(4.3) 



By (4.1), (4.3), (3.17a,b), 
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(4.4) 



By (4.1) and (4.2) we see that the Symanzik term is a surface term, which diverges as the 
cutoff is removed. However, it is there and, as remarked by [IH], cannot be disentangled from 
the finite part (4.4). 

However, we have a surprise: the same value of £f ((3.17b) or (4.4)) is obtained by Hadamard 
regularization of the divergent integral of E^"'^{z, d; D) given by (2.12a) - the resuh of Iqft - as 
we now show. Let P{.) denote the Hadamard regularization (or "partie finie") of a given integral 
(see ,pp. 38-43). Then 

Theorem 4.1 



where £/ is given by (3.17b). 
Proof. 



p(^J^E'""'{z,d;D)dz^ = Ef 



(4.5) 



dz- 



1 1 



hm 



4(27r)2 z^ J 4(27r)2 .^o 



P / dz 



dz 



4(27r)2(d-z)^ 



3e3 



1 1 

12(27r)2d3 



1 1 

12(27r)2d3 



(4.6a) 
(4.6b) 



which, inserted into (2.12), yields, after a simple calculation, 

which equals Sf by (3.17b). 

Remark 4.1 : (2.12a) shows that the divergence of the energy per unit area remaining after 
renormalization (1.5) - i.e., after performing the cancelation between the first integral in (2.3) 
with the m = term in Ii is independent of the ultraviolet cutoff, being only due to the sharp 
nature of the surface, i.e., the use of Dirichlet (or Neumann) b.c. on quantum fields. These are 
the "unnatural acts" referred to by R. Jaffe in In the Casimir problem this divergence is 
removed by Hadamard regularization, as proved in Theorem 4.1, yielding the correct energy per 
unit area. Such was also the finding of Elizalde in [21] : he showed that Hadamard regularization 
yielded, in models treated in [18], the same result for the finite part of the energy density. As 
seen in (4. 5)- (4. 7), this procedure identifies precisely the singularities associated to the infinitely 
thin (sharp) boundary. 



5. Conclusion 



20 



We have presented in this paper a rigorous method to study the Casimir energy density, 
and apphed it to the massless scalar field with Dirichlet b.c. on parallel plates. The massive 
case, as well as other geometries (like cube, sphere) may in principle be treated by this method, 
and independence of b.c. has been verified explicitly in each case (although no general theorem 
can be invoked). The finite part ej may be derived in two equivalent ways: 1) with the use of 
(ultraviolet regularization-dependent) Symanzik surface-renormalization counterterms, and 2) 
with no ultraviolet regularization, but applying Hadamard regularization to the energy density 
at the singular surfaces . In this connection, we remark that (1.6b) becomes a distribution after 
subtraction of a surface term which may be guessed from (4.1) and (4.2). Let PF denote the 
Hadamard pseudofunction [25](pg.41). By (2.12a), (2.13) and the distributional formula [25] (II 
1.28) 

i^PF(l)=PF(l) 

we obtain the final renormalized one point function 

coiz) = PF[E{z - ■)] (5.1) 

where E{z,d) is defined by 

-^H(z, d) = E'""'{z, d;D) iovO<z<d (5.2) 

and the dot in (5.1) stands for any point at any of the singular surfaces z = and z = d . 
Interestingly, this one point function (5.1), which is a constant for periodic b.c, has a singularity 
structure, with respect to the singular surfaces, of the same nature as the one displayed by the 
two-point function of a locally covariant quantum field theory [31j, i.e., of the form (see, e.g., 
1, and [30] 

1 rAi/2 1 

+ w (5.3) 



^2{xi,X2) = ^ PF 



Ai/2 

+ V Ina 



a 

where A^/^, v and w are certain smooth functions, and a = a{xi,X2) is the square of the 
geodesic distance between the points xi and X2- in (5.2) one has the square of the distance (here 
unambiguously defined) between the point z and any one of the singular surfaces. For the two- 
point function the Hadamard form has a fundamental significance |31] . and it seems worthwhile 
to inquire on a possibly deeper meaning of this structure in the present context. As remarked by 
de Witt [37], even the Casimir effect is primarily a problem of Riemannian manifold structure, 
and, in this sense, locally covariant quantum field theory might place some constraints on the 
problem. 

In conclusion, we have analysed a more realistic model - the massless scalar field enclosed by 
parallel plates - in the framework of Iqft introduced by B.Kay [9j. Two fundamental methods 
of asymptotic analysis play a dual role in our treatment. The Poisson summation formula |20j 
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yields a cutoff independent energy density, which is integrable and homogeneous in the case 
of periodic b.c, but is not integrable in the region between the plates in the Dirichlet case 
(also true for Neumann or mixed b.c). The Euler-MacLaurin summation formula shows the 
existence, in the asymptotic series for the energy per unit area, of a regularization dependent 
surface term, which diverges as the cutoff is removed, in the case of Dirichlet b.c: this term is 
absent if periodic b.c. are used (remark 3.6). 

We now summarize the final message of the present paper. In addition to local observables 
such as the energy density (2.12a,b), we propose to regard intensive variables such as the energy 
per unit area (2.18) as fundamental observables. This is motivated by the fact that in those cases 
in which the effect is measured, the observed quantity, the pressure, originates from (2.18) by the 
thermodynamic formula. In applications to Cosmology, for example, one is taken to calculate 
the Casimir energy of a "small" closed region (e.g., a cube) and the pressure by thermodynamics. 

Adopting this view, Iqft clearly selects periodic b.c- in this case ep is finite by (2.19b)- 
rejecting Dirichlet ( Neumann and mixed) b.c, in which case ed diverges by (2.19a). Another 
way to see this is by the emergence of a RDST in the latter, but not in the former, as shown in 
section 3. The emergence of the RDST may be seen as a rigorous proof, for the present model, 
of R.Jaffe's assertion |18j that the renormalized energy is, in general, material dependent, and 
diverges if a certain parameter controlling the "sharpness" of the surface (A in [18]) tends to the 
value corresponding to an infinitely thin or sharp surface (oo in [IB]). The RDST can certainly 
be made finite, in the limit where the cutoff is removed, by adopting methods such as those in 
[32], which allow a macroscopic description of a "soft" surface. The latter involves a nonlocal 
term, as in classical electrodynamics with dielectric materials. A solution of this problem by Iqft 
would probably require a detailed microscopic description of the surface (see, however, [17]). 

It may be asked whether the rejection of Dirichlet ( Neumann, mixed) b.c. by Iqft is phys- 
ically plausible, i.e., whether these b.c. are imposed by physics. The only known physical 
case is the e.m. theory with conductor surfaces (idealized as infinitely thin), where Dirichlet 
or Neumann b.c. impose themselves by reference to the classical limit, which occurs provided 
the particles constituting the conductor act collectively in an essentially classical manner [16], 
but in that case there is a miraculous cancellation ([4], [37]) and no problem arises! The scalar 
field counterpart to this assertion occurs in the conformally invariant scalar field, as mentioned 
in the introduction, in which case the nonintegrable terms in (2.12a) disappear, or, in different 
terms, the RDST in section 3 is absent, as happens in the ordinary scalar quantum field with 
periodic b.c. However, as argued in the introduction, the condition of tracelessness, essential 
for the vanishing of the surface terms, is not met by the (still open) applications of the Casimir 
effect in Cosmology. 

The surface, in Cosmology, is not made out of special materials, but rather used when 
attempting to restrict the field to some region of space-time which is not reduced to a point. 
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which is obviously required when deahng with quantum fields. The quantity analogous to e is 
the average of the energy density over some closed region. This may lead to a partial justification 
of the distinguished role of the periodic b.c. : they select just what we conjectured to be the 
"field component" of the vacuum energy ( equal to £f in (3.17b) up to a factor) which, moreover, 
may be positive in the massive case [36], and yields, through the rcnormalization prescription 
(1.5), a homogeneous energy density (2.19b) of the type which is required in general relativity. 
Note that, in the case of Dirichlet b.c, even the finite part of the energy density (5.1) is not 
homogeneous. A basic open problem is the dependence on the region's dimension (d^^ in 
(2.19b), for a cube of side L [6]): taking L of order of the "radius of the Universe" leads to 
an absurdly small value of the dark energy density [12]. 

Acknowledgements: We are greatly indebted to the three referees, most particularly refereel, 
whose comments considerably clarified the conceptual basis of this paper. 
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